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Abstract 

A rotational modulator (RM) gamma-ray imager is capable of obtaining significantly better angular resolution than the fundamental 
geometric resolution defined by the ratio of detector diameter to mask-detector separation. An RM imager consisting of a single grid 
of absorbing slats rotating ahead of an array of a small number of position-insensitive detectors has the advantage of fewer detector 
elements (i.e., detector plane pixels) than required by a coded aperture imaging system with comparable angular resolution. The 
RM therefore offers the possibility of a major reduction in instrument complexity, cost, and power. A novel image reconstruction 
technique makes it possible to deconvolve the raw images, remove sidelobes, reduce the effects of noise, and provide resolving 
power a factor of 6 - 8 times better than the geometric resolution. A 19-channel prototype RM developed in our laboratory at 
Louisiana State University features 13.8° full-angle field of view, 1.9° geometric angular resolution, and the capability of resolving 
sources to within 35' separation. We describe the technique, demonstrate the measured performance of the prototype instrument, 
and describe the prospects for applying the technique to either a high- sensitivity standoff gamma-ray imaging detector or a satellite- 
or balloon-borne gamma-ray astronomy telescope. 
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1. Introduction 

^ ■ Imaging hard x-ray and gamma-ray photons is not possi- 
C$ ble using traditional focusing techniques. Instead, an indirect 
l ™ — 1 technique must be used, such as spatial/temporal modulation or 
Compton scattering. A single-grid Rotational Modulator (RM) 
^ is one such instrument, which combines the advantages of two 
[more-commonly used imagers: the spatially-modulating coded 
aperture and the temporally-modulating double-grid Rotating 
^} Modulation Collimator (RMC). 

• ■ A coded aperture (Fig. [U left) (UBS] consists of a mask of 
opaque and transparent pixels (typically equal numbers of each, 
resulting in a 50% "throughput") arranged in a specific pattern. 
The source distribution in the object scene casts a shadow of the 
I* mask onto a position-sensitive detection plane. Angular resolu- 
tion of the instrument is given by the ratio of mask pixel size 
to mask-detection plane separation, and can be made very good 
in practice. However, the size of the detection plane elements 
must be small enough to allow for adequate measurement of the 
shadow (less than or equal to half the mask pixel size in order to 
satisfy the Nyquist criterion). For a high- sensitivity, large-area 
instrument, this requirement results in a high number of readout 
channels, increasing cost and complexity. 

A Rotational Modulation Collimator (Fig. [Q middle) HE] 
is one of a class of temporal modulators. It consists of two con- 
centric offset grids of slats above a single position-insensitive 
detector. As the grids rotate in tandem, photons from off-axis 
sources will be modulated before detection. The detector mea- 
sures a time history of counts, unique to the source distribution 
in the object scene. The angular resolution of the RMC, like 
the coded aperture, is defined by the width of the slat spacing 
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Figure 1: Basic structure of three multiplexing imagers. All instruments consist 
of a detection plane situated below a suspended mask which spatially or tempo- 
rally modulates incident photons. The RM shares its rotating grid in common 
with the RMC, and its single-grid and multiple detector elements with the coded 
aperture. 



divided by the mask-detector separation. Since the slat spac- 
ing can easily be made small, the RMC is capable of excellent 
angular resolution, despite having only a single readout chan- 
nel. The secondary grid, however, reduces the sensitivity of the 
instrument (25% throughput) and increases the overall weight. 

A Rotational Modulator (Fig. [T] right) (SLZD is another tem- 
poral modulation imager, but consists of only a single grid of 
opaque slats suspended above a small array of detectors. For 
a large-area gamma-ray detector, where the grid must be thick 
and correspondingly massive, the use of one rather than two 
grids can potentially be a significant advantage. The slat widths, 
slat spacing, and detector diameters are equal. The grid is ro- 
tated, periodically blocking photons incident on the instrument. 
This time-dependent shadow is recorded by the detectors as 
a history of counts during an entire exposure and then folded 
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Figure 2: (Top) Deterministic count profiles for a single point source viewed 
from three detector positions. The horizontal scale is grid rotation phase, while 
the vertical scale is unit source intensity. (Bottom) The point-spread function 
of the prototype LaBRAT RM imager described in the text in sec. 0below, with 
axes in degrees. 



modulo the grid rotational period to produce a single count pro- 
file for each detector. With temporal modulation and modest 
detection plane spatial resolution, the RM shares the advantages 
of both the coded aperture and the RMC. It maintains high sen- 
sitivity with a 50% throughput, low weight, and has a relatively 
simple readout system. Again, the intrinsic geometric angu- 
lar resolution is defined by the ratio of slat spacing to mask- 
detector separation. Since this spacing must be equal to the 
detector diameters, this geometric resolution is effectively lim- 
ited. Thus, for the advantages of the RM to fully outweigh those 
of other multiplexing instruments, the reconstruction technique 
must be capable of resolving images beyond the geometric res- 
olution, i.e. it must achieve "super-resolution." 

In the sections below, we briefly describe the iterative "Noise 
Compensating Algebraic Reconstruction" (NCAR) algorithm 
that has been developed to deconvolve RM images, the pro- 
totype RM imager that has been constructed in our laboratory 
at Louisiana State University (LSU), and the measurement re- 
sults that demonstrate the ability of the RM to image sources 
in the presence of background and to provide super-resolution 
several times better than the geometric resolution. We conclude 



by describing the design of an RM telescope suitable for high- 
sensitivity stand-off gamma-ray measurements or gamma-ray 
astronomy observations on a high altitude balloon platform. 



2. Image Reconstruction 

The RM grid rotates at a fixed frequency oj. An object scene 
S (n) produces a time history of counts in each detector d, 

O d (t) = YjP d (n,t)S(n) + N d (t). (1) 

n 

The instrument response Pd(n,t) is a collection of pre- 
calculated count profiles for all possible source locations n 
within the object scene. Fig. El (top) shows count profile exam- 
ples of the same source viewed by detectors at three different 
positions. 

To reconstruct the object scene, a cross-correlation image 
is first generated by the summed convolution of the measured 
count profiles from each detector with the instrument response: 

c ^ = Z Z Pd(n > t)0d(t) = Z w ' n)S (n/) + N(n) > (2) 

d t n' 

where 

P\n\ n) = Yj Z Pdin '> t)Pd ^ °' (3) 

d t 

The cross-correlation image for a point source reveals the PSF 
to have a ring pattern centered on the source location (Fig. [21 
bottom). Algebraic solutions attempt to solve directly for the 
object scene by deconvolution of the cross-correlation with the 
PSF. Direct Demodulation (DDM) 18D is one such technique, 
which solves the system of equations iteratively while enforcing 
a positivity condition as a physical constraint on the reconstruc- 
tion. DDM has been previously shown by its authors to provide 
super-resolution for RMCs |31, and we have shown the same 
ability with RMs in simulation where little or no noise exists 
fioh . Unlike a statistical technique, however, any algebraic re- 
construction method inherently requires that the reconstructed 
image agrees exactly with the data. Since the data include sta- 
tistical uncertainty from source and background measurements, 
this technique may lead to spurious source reconstructions and 
inaccurate true source locations when the measured signal-to- 
noise ratio (SNR) is low. 

NCAR is an iterative method devised to compensate for the 
presence of noise in the data BllN . The algorithm is based on 
the DDM technique but includes compensation for noisy data to 
suppress fluctuations. A modified Gauss-Seidel method is used 
to solve for the object scene 1 12]. A randomized background 
layer, derived from Poisson uncertainty in the data, is then 
added to the cross-correlation image at each iteration. Once the 
reconstruction converges to the data as well as possible given 
the noise limitations, subsequent reconstructions are combined 
in a running average; spurious random peaks are suppressed 
while true sources remain. 
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Figure 3: (Left) Photograph of the LaBRAT instrument. (Right) Closeup of 
LaBRAT's detection plane. 



3. Prototype RM Imager 



The Lanthanum Bromide-based Rotating Aperture Telescope 
(LaBRAT, Fig. [3} is a laboratory prototype RM developed at 
LSU fl3ll . It was designed to verify the feasibility of the RM 
technique and the reconstruction algorithm, and as a pathfinder 
experiment to an eventual high-altitude balloon mission for 
gamma-ray astronomy in the 30 - 700 keV energy range. While 
the RMC has seen limited use in gamma-ray imaging (e.g. 
[03]), the authors are aware of only one other single-grid RM- 
type instrument; the WATCH experiment [[TsL Eg], however, 
uses a unique and different system of detection than LaBRAT. 

The LaBRAT mask is a 2' diameter grid of eight 1.5" x 24" 
x 0.75" thick lead slats, spaced 1.5" apart. The slats are sup- 
ported by an aluminum structure, and sandwiched between two 
sections of PVC pipe. The grid sits atop four roller bearings to 
allow for free rotation. A servo motor is attached by a belt and 
drives the mask to rotate at approximately 10 rpm. Four pho- 
tointerrupters are mounted around the grid in order to orient the 
mask and allow for subsequent folding of the count profiles. 

The detection plane is composed of nineteen Cerium-doped 
Lanthanum Bromide (Ce:LaBr3) BrilLanCe 380 scintillators 
manufactured by Saint-Gobain [17]. LaBr3 provides excel- 
lent energy resolution compared to other inorganic scintillators 
(< 3% FWHM at 662 keV) due primarily to its high light yield. 
Each detector is hermetically sealed in an aluminum housing 
with a glass window on one end; this end is optically coupled 
to a 1.5" Electron Tubes 9102B photomultiplier tube. The de- 
tection plane is offset ~1.2 m from the mask, defining a 13.8° 
full-angle field of view and 1.9° geometric angular resolution. 

The data acquisition system was custom designed and built 
for LaBRAT. Each PMT anode signal is fed into a separate 
channel on the data acquisition front end module (FEM), where 
the signal is split; one of the signals goes to the pre-amplifier 
circuit and is eventually digitized, while the other is used for 
triggering. Currently, the system is capable of reading 32 chan- 
nels, but has been designed with a larger system in mind, suit- 
able for a balloon-flight instrument, and so is easily scaled. 
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Measured Count Profile (Outer Ring Detector) 




Figure 4: Normalized measured count profiles for a middle and outer-ring 
detector for one of the observations shown in Figure [5] 

4. Results 

Monte Carlo simulation results have been reported in fllll . 
Here, two 50 yuCi 137 Cs sources are used to demonstrate the 
measured resolving power of LaBRAT. The sources are placed 
-10 m from the detection plane. Since the RM was designed 
primarily as a far-field imager, the reconstruction algorithm as- 
sumes photons from a source are mutually parallel when inci- 
dent on the instrument. Parallax can degrade the image, and so 
this distance provides a suitable tradeoff between minimal par- 
allax and higher measured source rate. With a 5 Hz per source 
count rate and 30 Hz background in the selected energy range, 
exposure times of 5 - 10 hours are used for these measurements. 

The sources are placed at varying angular separations and 
imaged (Fig. 0. The 0.75" thick mask only attenuates ~ 90% 
of the 662 keV photons observed, and so an absorption term 
is included in the calculation of the count rate profile Pdif). 
At 2°22 / , the cross-correlation image shows the two sources, 
and the reconstruction successfully removes the sidelobes. At 
1°46', almost equal to the geometric resolution of the instru- 
ment, the cross-correlation image begins to deteriorate. At even 
smaller angular separations (1°10' and 35'), only a single cen- 
tral peak is visible in the cross-correlations. In all four cases, 
however, the final reconstructions derived from the NCAR al- 
gorithm are able to deconvolve the data and fully resolve the 
two sources down to a separation of at least 35', successfully 
achieving super-resolution of a factor greater than 3. As addi- 
tional improvements are made to the prototype, resolving power 
of 15 - 20' is expected, as has been shown in simulation 11 Oil . 

5. Discussion and Conclusions 

Coupled with an appropriate reconstruction technique, it ap- 
pears that a rotational modulator is a suitable alternative to a 
coded aperture. It maintains the same throughput, and is ca- 
pable of achieving comparable angular resolution with a far 
simpler and cheaper detection plane. We have built and tested 
a laboratory prototype, LaBRAT, to demonstrate this concept, 
and have shown the ability to achieve angular resolution ~ 3x 
better than the intrinsic geometric angular resolution defined by 
the instrument. 
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Figure 5: Contour plot of cross-correlation image (upper row), contour plot 
of final NCAR reconstruction (middle row), and surface mesh of final recon- 
struction (bottom row) for two sources at four different angular separations, as 
imaged by LaBRAT. Note, in particular, the cross-correlation images for 1°10' 
and 35' separation, which are not capable of resolving the two sources inde- 
pendently; the final reconstructions using the NCAR algorithm (middle and 
lower images), however, are capable of achieving this resolving power. Axes in 
degrees. 



Calibrations of the newly-built data acquisition system are 
still ongoing. Based on earlier simulations and results with a 
previous version of the data acquisition system [ 18], we expect 
to be able to separate sources at 15 - 20'. Further improvement 
is expected by dithering the telescope around the source posi- 
tion. (For a balloon-borne version of the instrument used for 
gamma-ray astronomy, as described below, the motion of the 
balloon platform will provide a natural dithering motion.) A 
detailed discussion of the effects of dithering, the applications 
to moving sources and moving detectors, and angular resolu- 
tion results as a function of signal-to-background rates will be 
published in a separate report. 

The High Altitude Rotational Modulator for Energetic radia- 
tion Imaging (HARMEnI) gamma-ray astronomy telescope de- 
sign is based on the LaBRAT prototype results and the NCAR 
imaging algorithm. HARMEnI will serve as a pathfinder mis- 
sion for a satellite-borne instrument to survey and monitor 
black-hole candidates, observing at energies between 30 and 
700 keV. Relative to LaBRAT, the initial HARMEnI telescope 
will have increased sensitivity and field of view. Ten lead bars 
comprise the mask, and the detection plane is a tightly-packed 
arrangement of 37 LaBr3 crystals - the same as those used with 
LaBRAT. HARMEnI provides a field of view of 20° and geo- 
metric angular resolution of 1.9°. Since HARMEnI' s geometric 
resolution is the same as that of the LaBRAT prototype, an ul- 
timate resolving power of 15 - 20' is anticipated. 

Background rates at balloon float altitude have been mea- 
sured with a similar Ce:LaBr3 crystal flown on a long duration 
(17 day) Antarctic balloon flight [19]. The main contributors 
to the background are the atmospheric and cosmic diffuse flux 
and the internal background due to 138 La decays. The expected 
background rates for HARMEnI are estimated using the mea- 
sured data from the Antarctic flight and assuming a 5 mm thick 
passive graded lead-tin-copper shield extends from the detec- 
tion plane to the mask and covers the back of the detectors. 
Based on the expected rates, the observation time necessary for 
the Crab Nebula will be significantly less than one hour. 
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